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High-resolution scanning Hall probe microscopy has been used to image vortex configurations in 
very large periodic arrays of artificial pinning sites. Strong matching effects are seen at fields where 
either one or two vortices can sit at a site; with three vortices per site, however, no clear matching 
is observed. Matching effects have been also been observed at several fractional multiples of the 
matching field, including 1/5, 1/4, 1/3, 1/2, and 3/4. These fractional values are characterized by 
striking domain structure and grain boundaries. 
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The behavior of superconducting vortices in the pres- 
ence of a periodic array of holes reveals rich and unex- 
pected static and dynamic phenomena. Transport 0-|^] 
and magnetization studies show distinctive features 
at matching fields where the vortex structure is commen- 
surate with the hole array. Commensurability can arise 
from vortex configurations [^-^) in which each hole con- 
tains an integer number of flux quanta, or from intersti- 
tial structures where vortices occupy the regions between 
flux-saturated holes. Interesting configurations also arise 
at fractional matching fields, where the occupation num- 
ber of each hole or the number of interstitial vortices 
forms a superstructure locked to the basic hole array . 

Although magnetization and transport studies have 
elucidated much of the basic phenomenology of the ar- 
ray/vortex system, they can measure only its global prop- 
erties and cannot deduce details of the local configura- 
tional vortex state. More recently Lorentz microscopy 
has been used to image vortex configurations || , but only 
in relatively small arrays consisting of 13 x 13 = 169 pin- 
ning sites. In this paper we present large-area scanning 
Hall probe microscope (SHM) images of vortex config- 
urations in arrays containing sa 10 6 holes. Our images 
span some 5000 holes, and so yield important informa- 
tion on the large-scale structure of the vortex configura- 
tions. These studies reveal striking multi-quantum and 
interstitial vortex patterns in a square-periodic hole ar- 
ray. At fractional matching fields we resolve distinctive 
domains of phase-slip related vortex superstructures that 
are separated by domain walls with characteristic inter- 
nal structure. 

The SHM uses a large-range scanner yielding « 130 
x 130 /im images. The sample investigated was a 100- 
nm-thick niobium film with 0.3-/im-diameter holes on 
a square lattice. The lattice constant a of the film 
was 1.870 /mi, yielding an expected first matching field 



H m = &o/a 2 = 5.913 G. The array was produced using a 
lithographic technique based on the interference of light 
, which yields very large and uniform arrays (here « 2 
mm x 3 mm). Near T c — 8.37 K, the array exhibits clear 
cusps in its magnetization at H m and 2H m , but not at 
any higher integral multiples of H m . 

SHM images near several matching fields are shown 
in Fig. 1. For each image, the samples were field-cooled 
from above T c to a base temperature below 3 K where the 
images were taken. The top row of Fig. 1 shows the pro- 
gression in the vortex configurations as the applied field 
H is varied around H = H/H m — 0. At H w we see 
several isolated vortices (dark spots) . The actual vortex 
diameter is much smaller than these spots, whose size is 
determined by the Hall probe resolution. A somewhat 
darker spot, just to the right and below center, appears 
in all scans and is presumably associated with a physical 
defect in the hole lattice, which allows extra flux to sit at 
that point. As the field is increased from zero, vortices 
begin to enter the sample in larger numbers, as seen in 
the H = 0.075 and 0.15 images. A similar progression 
is observed for negative field increments away from zero, 
as shown at H = —0.075 and —0.15. Here the vortices 
appear as light spots, indicating that they arc oppositely 
directed to those seen in positive fields. 

As we continue to increase the field towards what we 
will be able to identify as the first matching field H m 
(H = 1), we begin to see a very remarkable progression 
(Fig. 1, second row). For instance, at an applied field 
of H = 0.85 = 1 — 0.15, the image looks statistically 
identical to the image at H = —0.15! And, remarkably, 
the image at H = 1 is essentially indistinguishable from 
that at H = 0. It is important to note here that in each 
image the average value of the field has been subtracted 
out. This is why the H = 1 image at about 6 G and the 
H = image appear to have about the same overall gray 
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level. We can then understand the appearance of the 
H = 1 scan as follows. The held increment from an ex- 
actly H = image (no spots) to the next no-spot image 
is 5.929 G, very close to the hrst matching held of 5.913 G 
deduced from the lattice geometry. Thus we identify this 
later spot-free image at 5.929 G with the hrst matching 
held or H = 1. The smooth gray background near H « 1 
can then be explained in the following way. Exactly at 
H = 1, each hole contains exactly one vortex, leading 
to a dense and uniform configuration of vortices. Evi- 
dently the spatial resolution of our Hall probe is not good 
enough to image individual vortices when they are this 
close together and perfectly ordered fl(i|| , so this regular 
array of vortices appears as a smooth gray background. 
The few black spots in the H w 1 image are thus (easily 
visible) extra vortices. 

By the same token the white spots in the images at 
H = 0.85 and 0.925 are not antivortices, as they were 
in the —0.15 and —0.075 images. Instead, they are va- 
cancies in a smooth background of filled holes. Since 
vacancies — holes with no vortex — have a lower (whiter) 
held than the broad gray areas which are filled with vor- 
tices, they appear as white dots. 

The third row of Fig. 1 shows the vortex configura- 
tions as the held is further increased through the second 
matching held (H — 2). A very similar progression is 
observed. At H w 2, there is again a smooth background 
populated by a few black spots. Again, we interpret this 
background as the highly ordered state with now two vor- 
tices sitting in each hole; black spots for H > 2 are again 
extra vortices, and white spots for H < 2 now represent 
a hole with only one vortex. At the third matching held 
(H = 3), however, the appearance is radically different 
(Fig. 1, bottom right). Instead of the smooth gray back- 
ground we would expect if three vortices sat in each hole, 
we see a rather muddled image with no discernible struc- 
ture. We believe this different appearance results from 
the sudden presence of interstitial vortices above H = 2. 

To probe this issue more directly, we have taken 
closeup scans just below and above H = 2. Figure 2(a) 
shows a 25 /j,m x 28 /jm scan at a field H = 2 — 0.084, 
and Fig. 2(b) shows ~H = 2 + 0.084. Also shown is the 
least-squares-fit positions of the holes as determined from 
an image taken at H = 1/2 where the hole positions are 
unambiguous [see, e.g., Fig. 3(b)]. We also have used 
an absolute position sensor |pd[| mounted on the scanner 
head to compensate for any possible drifts between im- 
ages. In Fig. 2(a) each white spot (i.e., a "missing" vor- 
tex) clearly sits directly on a hole, indicating as well that 
all vortices also sit on the holes. However, in Fig. 2(b) the 
black spots, which are now extra vortices, sit variously 
on holes or on interstitial sites. As is well known, each 
hole in a him can hold up to a saturation number n s of 
vortices ]l2[ ]; with our material parameters we estimate 
that n s is near the boundary between 2 and 3. With 
the sputtering/lift-off lithography used here M there will 



be inevitable fluctuations in the hole diameter, perhaps 
yielding n s = 2 for some holes and n s = 3 for others. 
Vortices nucleating near n s = 2 holes would end up as 
interstitials, while those near n s = 3 sites could sit in a 
hole. This random-appearing admixture of vortices on 
holes and in interstitial sites would, by H = 3, lead to 
the very disordered configuration observed at that field 
(Fig. 1). 

It is also possible to have matching effects at non- 
integer multiples of H m . Indeed, distinct features at 
H = 1/4, 1/2, and possibly 1/5 (or 3/16), 1/8, and 1/16 
have been observed in magnetization studies ||,[n|], and 
Harada et al. f| have imaged non-integer matching at 
1/4, 1/2, 3/2, and 5/2. In our imaging experiments, we 
also observe matching effects at several of these fractional 
fields, as well as some not previously observed. We find 
that vortex configurations at fractional matching fields 
are characterized by striking domain structure and asso- 
ciated grain boundaries. 

By far the strongest submatching effects occur at H = 
1/2. Figure 3(a) shows that the configuration consists of 
large areas of well-matched vortices separated by curious 
stripe-like boundaries. The closeup in Fig. 3(b) reveals 
that the smooth regions consist of vortices occupying ev- 
ery other hole in a checkerboard-like fashion. Regions 
of different "polarity" — with the vortices occupying ei- 
ther the "red" or "black" squares of the checkerboard — 
are separated by striped grain boundaries [B. Fig. 3(c) 
schematically shows how a polarity shift results in rows 
or columns of alternating pairs of vortices and vacancies, 
which appear in the images as the striking boundary fea- 
tures. At the junctions between north-south and east- 
west boundaries there is always a bright or dark spot, 
whose origin is again clear from Fig. 3(c). From con- 
siderations of the boundaries, we have mapped out the 
domains of differing polarity [inset, Fig. 3(a)]. One cu- 
rious feature of these grain boundaries is that they run 
predominately north-south, indicating a breaking of the 
square symmetry of the lattice in these two directions. 

Weaker but still striking matching effects occur at 
H = 1/4 as well as its complimentary field of 
~H = 1 - 1/4 = 3/4 (Fig. 4). We note immediately that 
the matching here is much less distinctive than at 1/2, 
1, or 2. Nonetheless, there do exist some relatively large 
ordered domains. Zooming in on the boundary of such 
a patch [outlined in Figs. 4(a) and 4(b)] allows us to ex- 
plore the vortex configuration in detail [Figs. 4(c) and 
4(d)]. Here we can clearly see areas with well-ordered 
vortices interspersed with phase slips (diagonal light and 
dark segments) and other, more complex features. The 
vortex conhguration in the ordered regions is shown to 
scale schematically in Fig. 4. Notice that the configura- 
tion may be viewed in terms of alternating empty and 
half- filled rows ||[l6| ; this again indicates the presence of 
a symmetry-breaking field which selects this particular 
grain orientation. One possible source of this asymme- 
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try could be a small difference in the interhole spacing 
a in the horizontal and vertical directions. However, no 
such asymmetry is visible at the w 0.1% accuracy of our 
optical diffraction measurements of a. 

Besides the 1/2, 1/4, and 3/4 states already discussed, 
we have found several other relatively well-ordered states 
at fractional fillings < H < 1. Figure 5(a) shows the 
configuration at H = 1/5. Although there are no large 
domains, in much of the image we can observe a gen- 
eral tendency towards forming a structure which con- 
sists of a square lattice of vortices tilted by an angle 
arctan(l/2) w 26.6° p|,|i"6||; this structure is shown in 
detail in the inset. There is also striking structure at 
H = 1/3 [Fig. 5(b)]. Here, the configurations consist of 
slope-one diagonal stripes of (unresolved) vortices sepa- 
rated by two empty diagonal stripes. Such a configura- 
tion has been predicted by Watson |l6| for the general 
case of repulsive particles on a square lattice, and by 
Teitel and Jayaprakash jl7| for the related case of vor- 
tices in a superconducting wire network |f[8| . Again, do- 
mains with this precise structure are small, and there 
are a large number of defects present. There is very 
weak ordering present at H = 2/5 [Fig. 5(c)], where 
Teitel and Jayaprakash predict states consisting of 
two empty diagonals and then a full-empty-full diago- 
nal sequence. This later sequence has the appearance of 
a checkerboard state, and is barely visible between the 
bright (empty) diagonal stripes of the inset to Fig. 5(c). 
Finally, we show an image taken "at" the irrational field 
H = (V5 — l)/2 ps 0.618. No apparent order is present 
in this case. 

Vortex configurations in square-periodic hole arrays 
thus reveal remarkably complex patterns reflecting the 
interplay between the pinning energy of the hole array 
and the interaction energy of the multi-quantum and in- 
terstitial flux structures. Our images show clear multi- 
quantum occupation of the holes up to saturation, fol- 
lowed by the appearance of interstitial Abrikosov vortices 
at higher fields. At fractional filling factors, domains 
of well-ordered vortices are separated by a network of 
domain boundaries displaying remarkably simple order 
for half filling and increasing complexity at other fillings. 
There appears to be a nontrivial set of ordered fractional 
states, a complete inventory of which awaits a general 
theory correctly incorporating the several competing in- 
teractions present. Our Hall probe images thus provide 
direct visual evidence for a rich diversity of collective 
ground states arising from competing commensurability, 
order, and randomness in the vortex phase diagram of 
periodic hole arrays. 
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FIG. 1. SHM images obtained at applied magnetic fields H 
near the matching fields H = H/H m = 1, 2, and 3. The scans 
are 124 /im x 138 /im, contain some 5000 holes, and have a full 
scale of 0.73 G. Note the striking similarity between the vortex 
configurations near H — 1 and 2 and those near H = 0. The 
smooth backgrounds at H — 1 (2) consist of holes uniformly 
filled with 1 (2) vortices. At H — 3, however, the vortex 
configuration is highly disordered because of weakly-pinned 
interstitials competing with vortices in holes. The dark black 
spot is a local defect in the sample. 



FIG. 2. Vortex configurations (a) just below (H = 1.916) 
and (b) just above (H = 2.084) the second matching field. 
The small dark circles are the positions of the holes. Below 
H — 2 the vacancies (white spots) all sit directly on holes; 
thus all vortices must as well. Above H = 2, the extra vortices 
(dark spots) sit on both hole and interstitial sites. 



FIG. 3. (a) Vortex configuration at H = 1/2. There are 
large areas of vortices arranged in a checkerboard pattern; 
domains [outlined in inset to (a)] of opposite checkerboard po- 
larity are separated by stripe-like grain boundaries. A closeup 
of this configuration is shown in (b) and schematically in (c). 
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FIG. 4. Vortex configurations at H = 1/4 (a) and 3/4 (b). 
Selected regions are expanded in (c) and (d). The config- 
urations consist of relatively small well-ordered regions sur- 
rounded by complex disordered boundaries. Also shown are 
the inferred vortex arrangements in the ordered regions. Light 
circles represent empty holes, and dark circles occupied holes. 

FIG. 5. Vortex configurations at several fractional values 
of H. Also shown is the "irrational" field (f> = (V§ — l)/2. 
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